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Abstract 

Recent progress in the calculation of radiative correc- 
tions and in Monte Carlo event generation, relevant for a 
future e+e~ linear collider, is reviewed. 

INTRODUCTION 

Precision measurements at LEP, SLC, and the Tevatron 
rendered the last decade the era of high-precision physics. 
A future c + c~ linear collider (LC), such as TESLA [1], the 
NLC [2], or the GLC (former JLC) [3], does not only offer 
an even greater physics potential, but in turn represents a 
great challenge for theorists to understand phenomena at 
the experimentally achievable level of precision. 

For instance, returning again to the Z-boson resonance in 
the "GigaZ" mode of TESLA (where about 10 9 Z bosons 
can be produced within 50-100 days of running) allows 
for a repetition of the LEP1/SLC physics program with 
roughly an order of magnitude higher precision (see also 
Ref. [4]). Specifically, the uncertainty in the effective 
weak mixing angle could be reduced from 1.7 x 10~ 4 to 
1.3 x 10~ 5 . For a theoretical description of the Z reso- 
nance at this level of accuracy full two-loop calculations 
of the observables as well as the knowledge of leading 
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higher-order effects are clearly necessary. A scan over 
the W-pair production threshold could provide a sensitiv- 
ity to the W-boson mass of about 7MeV, which should 
be compared with the present error of 34 MeV, resulting 
from the kinematical reconstruction of W bosons. The 
present approach of approximating the radiative correc- 



tions to e + e 



WW — > 4fermions by an expansion 



about the double resonance is not applicable (or at least 
not reliable) in the threshold region where singly- or non- 
resonant contributions become important. The only solu- 
tion seems to be the full treatment of the complete four- 
fermion production processes at the one-loop level, includ- 
ing higher-order improvements. 

At energies exceeding the reach of LEP2, many new 
processes will be accessible, such as top-quark pair pro- 
duction, Higgs production (if the Higgs boson exists), or 
reactions with new-physics particles, as e.g. predicted by 
SUSY models. Most of these heavy particles are unsta- 
ble, so that their production eventually leads to many- 
particle final states. For example, the production of tt pairs 
or of a Higgs-boson with an intermediate or large mass 
(Mh > 2Mw) leads to six-fermion final states. To ex- 
ploit the potential of a LC, predictions for such reactions 
should be based on full transition matrix elements and im- 
proved by radiative corrections as much as possible. The 
higher level of accuracy at a future LC does, however, not 
only call for proper event generators for many-particle fi- 
nal states. "True" event generators, i.e. including parton 
showering and hadronization, have to be improved as well. 

In this brief article the main progress on precision cal- 
culations and event generators that has been achieved in 
the "Loopverein" and "Generators" working groups of the 
Extended ECFA/DESY Study since the appearance of the 
TESLA TDR [1] is reviewed. More studies on the physics 
potential of a LC in view of electroweak and strong inter- 
actions, top-quark physics, Higgs physics, and new physics 
searches, in particular supersymmetry, are summarized in 
Ref. [5]. 



HIGH-PRECISION OBSERVABLES AND 
MULTI-LOOP CALCULATIONS 

Precision calculations for n decay 

The precision measurement of the muon lifetime, or 
equivalently of the Fermi constant G M , sets an important 
constraint on the SM parameters, 



G,= 



V2M, 



2 „2 
Ww 



(1 + Ar), 



(1) 



where s 2 = 1 — = 1 — M^/M| and the quantity Ar 
comprises the radiative corrections to muon decay (apart 
from the photonic corrections in the Fermi model). In the 
past it has become common practice to implicitly solve this 
relation for the W-boson mass Mw> thus yielding a precise 
prediction for M\y that can be compared with the directly 
measured value. Recently the full prediction at the two- 
loop level has been completed. In detail, first the contribu- 
tions from closed fermion loops and from bosonic loops in- 
volving Higgs-boson exchange were calculated in Ref. [6] 
by making use of the FeynArts package [7] and the 
program TwoCalc [8] both written in Mathematica. 
The algebraic reduction leads to two-loop master integrals 
which are evaluated by semi-analytical methods. The full 
bosonic corrections have been calculated in Refs. [9] and 
[10]; Ref. [9] includes also a recalculation of the fermion- 
loop correction. In the former calculation the diagrams 
were generated with the C++ library DlAGEN (by Czakon) 
and evaluated using semi-analytical methods. In the latter 
case the graphs were generated with the package Diana 
[11] and evaluated by asymptotic expansions. The results 
of Refs. [6, 9, 10] are in good numerical agreement [12]. 1 
The two-loop fermionic corrections influence the Afw pre- 
diction at the level of ~ 50 MeV, where the bulk of this 
effect is due to universal, top-mass enhanced corrections to 
the p-parameter, which are proportional to m\ or m 2 . The 
non-universal two-loop fermionic corrections have an im- 
pact of up to 4 MeV, the two-loop bosonic corrections of 
only 1-2 MeV. 

The predictions at the two-loop level have been further 
improved by universal higher-order corrections to the p- 
parameter. The corrections of 0(G 2 m 4 a s ) and C(G^m^ ) 
have been calculated for arbitrary Mh in Ref. [13] (for 
other universal corrections to Ap and Ar see references 
therein) and were found to change Afw at the level of 
5 MeV and 0.5 MeV, respectively. The Feynman diagrams 
were generated using Qgraf [14] and asymptotically ex- 
panded with the program Exp [15]; the resulting massive 
three-loop tadpole integrals were evaluated with Matad 
[16]. 

Figure 1 compares the prediction for Mw, including the 
above-mentioned two-loop and leading three-loop effects, 
with the experimental value. Note that the shown paramet- 
ric uncertainty is much larger than the estimated theoretical 
uncertainty, which is about 3—4 MeV [17, 18]. Comparing 
this estimate with the aimed precision of 7 MeV in the A% 
determination at a future LC, the prediction of the W-boson 
mass from muon decay is in rather good shape. 

Precision observables on the Z resonance 

In order to describe the Z-boson resonance at LEP1 
within satisfactory precision it was possible to parametrize 
the cross section near the resonance in such a way [19] 
that a Born-like form with generalized "effective" cou- 

'The results of Ref. [6] for Afw have been corrected at the level of 
~ 1 MeV recently. 
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Figure 1 : Theoretical prediction of the W-boson mass with 
the parametric error from uncertainties in the top-quark 
mass and the running electromagnetic coupling, in com- 
parison with the experimental value MSZ P (plot shown in 
Ref. [17]) 



plings is convoluted with QED structure functions mod- 
eling initial-state radiation (ISR). From these effective Z- 
boson-fermion couplings so-called "pseudo-observables" 
were derived, such as various asymmetries, the hadronic Z- 
peak cross section, partial Z-decay widths, etc. Following 
the formal tree-level parametrization of the couplings, an 
"effective weak mixing angle", usually given as sin 2 #/, e ff, 
was derived for each fermion. Among these parameters 
the leptonic variable sin 2 (9i op . ff plays a particularly im- 
portant role, since it is measured with the high accuracy 
of 1.7 x 10~ 4 and is very sensitive to the Higgs-boson 
mass. The state-of-the-art in the precision calculations of 
the pseudo-observables, which is implemented in the pro- 
grams ZFITTER and TopazO (see Ref. [20] and references 
therein), did not change very much since the release of the 
TESLA TDR [1]. For instance, the estimated theoretical 
uncertainty in sin 2 6?i e p,eff is still ~ 6 x 10 -5 . A criti- 
cal overview about high-precision physics at the Z pole, 
in particular focusing on the theoretical uncertainties, can 
be found in Ref. [21] (see also Ref. [4]). 

Whether the pseudo-observable approach will also be 
sufficient for Z-boson physics at the high-luminosity GigaZ 
option remains to be investigated carefully. In any case, 
tremendous theoretical progress will be needed to match 
the aimed GigaZ precision on the theoretical side. For ex- 
ample, the expected experimental accuracy in sin 2 #i e p,eff 
is about 1.3 x 10~ 5 , i.e. about a factor 4 below the present 
theoretical uncertainty. A full control of observables at the 
two-loop level, improved by leading higher-order effects, 
seems to be indispensable. 



Recent results from the 2-loop frontier 

Although there are no complete next-to-next-to-leading 
(NNLO) predictions for 2^2 scattering reactions and 
1 — ► 3 decays (with one truly massive leg) available yet, 
enormous progress was reached in this direction in recent 
years. 

Complete virtual two-loop amplitudes for (massless) 
Bhabha scattering [22], light-by-light scattering [23], and 
e + e~ — ► 3 jets [24] have been worked out, using a large va- 
riety of special techniques, which have been summarized in 
Ref. [25]. A survey of similar results relevant for hadron- 
collider physics can also be found there. Apart from this 
two-loop progress on massless particle scattering, also a 
first step has been made towards massive Bhabha scatter- 
ing in Ref. [26]. 

Full NNLO calculations have to include real double- 
parton bremsstrahlung as well as interference contributions 
of one-parton bremsstrahlung and one-loop diagrams. The 
major complication in these parts concerns the proper ex- 
traction of the infrared (soft and collinear) singularities. 
The general form of multiple singular particle emission has 
been worked out in Ref. [27], which can serve as a basis 
for the extraction of the singularities. The actual separa- 
tion of the singularities can be performed either by apply- 
ing phase-space cuts ("slicing approach") or by subtract- 
ing an auxiliary cross section with the same singular struc- 
ture as the original integrand ("subtraction approach"). In 
Ref. [28] subtraction terms have been constructed for the 
leading colour contribution to c+e~ — > 2jets in NNLO. 
However, suitable general subtraction terms as well as their 
integrated counterparts that have to be added again are not 
yet available. 

Electroweak radiative corrections at high ener- 
gies 

Electroweak corrections far above the electroweak scale, 
e.g. in the TeV range, are dominated by soft and collinear 
gauge-boson exchange, leading to corrections of the form 
a N ln M (s/M^) with M < 2N. The leading terms (M = 
2N) are called Sudakov logarithms. At the one-loop (N = 
1) and two-loop (N = 2) level the leading and subleading 
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virtual and real corrections, these terms need not compen- 
sate in the electroweak SM for two reasons. The weak 
charges of quarks, leptons, and electroweak gauge bosons 
are open, not confined, i.e. there is (in contrast to QCD) no 
need to average or to sum over gauge multiplets in the ini- 
tial or final states of processes. Even for final states that are 
inclusive with respect to the weak charges Sudakov loga- 
rithms do not completely cancel owing to the definite weak 
charges in the initial state [30]. Moreover, the large W- 
and Z-boson masses make an experimental discrimination 
of real W- or Z-boson production possible, in contrast to 
unobservable soft-photon or gluon emission. 

In recent years several calculations of these high-energy 
logarithms have been carried out in the Sudakov regime, 
where all kinematical invariants (piPj) of different particle 
momenta pi, pj are much larger than all particle masses. 2 
A complete analysis of all leading and subleading loga- 
rithms at the one-loop level can be found in Ref. [31]. Dia- 
grammatic calculations of the leading two-loop Sudakov 
logarithms have been carried out in Refs. [29, 32]. Di- 
agrammatic results on the so-called "angular-dependent" 
subleading logarithms have been presented in Ref. [29]. 
All these explicit results are compatible with proposed 
resummations [33, 34] that are based on a symmetric 
SU(2)xU(l) theory at high energies matched with QED 
at the electroweak scale. In this ansatz, improved matrix 
elements M result from lowest-order matrix elements Mo 
upon dressing them with (operator-valued) exponentials, 



M ~ Mo ® exp (S ew ) <g> exp (S e 



(3) 



revealing that these corrections become significant in the 
high-energy phase of a future LC. In contrast to QED and 
QCD, where the Sudakov logarithms cancel in the sum of 



Explicit expressions for the electroweak and electromag- 
netic corrections 6 CW and 5 cm , which do not commute with 
each other, can, for instance, be found in Ref. [29]. For 
2^2 neutral-current processes of four massless fermions, 
even subsubleading logarithmic corrections have been de- 
rived and resummed [34] using an infrared evolution equa- 
tion that follows the pattern of QCD. 

In supersymmetric models the form of radiative correc- 
tions at high energies has also been worked out for a broad 
class of processes [35]. Based on one-loop results their ex- 
ponentiation has been proposed. 



Higher-order initial-state radiation 

Photon radiation off initial-state electrons and 
positrons leads to large radiative corrections of the form 
a N hi N (m 2 / s) . These logarithmic corrections are uni- 
versal and governed by the DGLAP evolution equations. 
The solution of these equations for the electron-photon 
system yields so-called structure functions, generically 
denoted by T{x) below, which can be used via convolution 
to improve hard scattering cross sections <r(p e +,p e -) by 



2 Note that this regime does not cover the case of forward scattering of 
particles, which is also of interest in several cases. 



photon emission effects, 



^(Pe+,Pe-) = / dx+T(x+) / da;_r(ar_) 
Jo Jo 

X &(X+P e +,X-P e -)- (4) 

While the soft-photon part of the structure functions 
(x — > 1) can be resummed, resulting in an exponential 
form, the contributions of hard photons have to be calcu- 
lated order by order in perturbation theory. In Ref. [36] the 
structure functions are summarized up to 0(a 3 ). Ref. [37] 
describes a new calculation of the (non-singlet) contribu- 
tions up to 0(a 5 ) and of the small-x terms [aln 2 (a;)] A ' 
to all orders (for previous calculations see papers cited in 
Ref. [37]). 

RADIATIVE CORRECTIONS TO 
2 3, 4, . . . PROCESSES 

W-pair production and four-fermion final states 

The theoretical treatment and the presently gained level 
in accuracy in the description of W-pair-mediated 4/ pro- 
duction were triggered by LEP2, as it is reviewed in 
Refs. [36, 38]. The W bosons are treated as resonances in 
the full 4/ processes, e + e~ — > 4/ (+ 7). Radiative correc- 
tions are split into universal and non-universal corrections. 
The former comprise leading-logarithmic corrections from 
ISR, higher-order corrections included by using appropri- 
ate effective couplings, and the Coulomb singularity. These 
corrections can be combined with the full lowest-order ma- 
trix elements easily. The remaining corrections are called 
non-universal, since they depend on the process under in- 
vestigation. For LEP2 accuracy, it was sufficient to in- 
clude these corrections by the leading term of an expansion 
about the two W poles, defining the so-called double-pole 
approximation (DPA). Different versions of such a DPA 
have been used in the literature [39, 40, 41]. Although 
several Monte Carlo programs exist that include universal 
corrections, only two event generators, YFS WW [40] and 
RacoonWW [41, 42], include non-universal corrections. 

In the DPA approach, the W-pair cross section can be 
predicted within ~ 0.5%(0.7%) in the energy range be- 
tween 180 GeV (170 GeV) and - 500 GeV, which was 
sufficient for the LEP2 accuracy of ~ 1% for energies 
170-209 GeV. At threshold (y/s < 170 GeV), the present 
state-of-the-art prediction results from an improved Born 
approximation based on leading universal corrections only, 
because the DPA is not reliable there, and thus possesses 
an intrinsic uncertainty of about 2%. In Figure 2 this un- 
certainty is compared with the sensitivity of the W-pair 
production cross section to the W-boson mass and some 
assumed experimental data points of a threshold scan. 
The figure demonstrates the necessary theoretical improve- 
ments. At energies beyond 500 GeV effects beyond 0(a), 
such as the above-mentioned Sudakov logarithms at higher 
orders, become important and have to be included in pre- 
dictions at per-cent accuracy. 
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Figure 2: Sensitivity of the W-pair production cross section 
to the W-boson mass and some assumed experimental data 
points, compared with the theoretical uncertainty of ~ 2% 
(taken from Ref. [1]) 
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At LEP2, the W-boson mass is determined by the re- 
construction of the W bosons from their decay products 
with a final accuracy of about 30McV. In Ref. [43] the 
theoretical uncertainty is estimated to be of the order of 
~ 5 MeV, which is comparable to the estimated [44] ac- 
curacy of ~ 10 MeV at a future LC. Theoretical improve- 
ments are, thus, desirable. 

The main sensitivity of all observables to anomalous 
couplings in the triple-gauge-boson vertices is provided by 
the W-pair production angle distribution. Figure 3 shows 
the impact of an anomalous coupling A 7 = ±0.03, the 
size of which is of the order of the LEP2 sensitivity, to- 
gether with the impact of non-universal corrections on the 
spectrum. The theoretical uncertainty in constraining the 
parameter combination A = A 7 = Az was estimated to 
be about 0.005 [46] for the LEP2 analysis. Since a future 
LC is more sensitive to anomalous couplings than LEP2 by 
more than an order of magnitude, a further reduction of the 
uncertainties by missing radiative corrections is necessary. 
However, a thorough estimate of the theoretical uncertainty 
in the determination of anomalous couplings at higher scat- 
tering energies (y/s > 500 GeV), where the experimental 
sensitivity to non-standard couplings increases, is not yet 
available. 

The above discussion illustrates the necessity of a full 
one-loop calculation for the e + e~ — > 4/ process and of 
further improvements by leading higher-order corrections. 

Single-W production 

The single-W production process e + e~ — > ev c W — * 
cv c + 2/ plays a particularly important role among the 
4/ production processes at high scattering energies. The 
process is predominantly initiated by 07* collision (see 
Figure 4) where the photon is radiated off the electron (or 
positron) by the Weizsacker- Williams mechanism, i.e. with 
a very small off-shellness qi. Consequently the cross sec- 
tion rises logarithmically with the scattering energy and is 
of the same size as the W-pair production cross section at 
about ^fs = 500 GeV; for higher energies single-W domi- 
nates over W-pair production. 

Theoretically the dominance of photon exchange at low 
q^ poses several complications. Technically, q^ — ► 
means that the electrons (or positrons) are produced in the 
forward direction and that the electron mass has to be taken 
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Figure 4: Generic diagram for the dominant contributions 
to single-W production, e + e~ — ► e~^ c W — > e~v e + 2/ 



into account in order to describe the cross section there. 
Moreover, the mere application of s-dependent leading- 
logarithmic structure functions does not describe the lead- 
ing photon-radiation effects properly, since ISR and final- 
state radiation (FSR) show sizeable interferences for for- 
ward scattering. Thus, the improvement of lowest-order 
calculations by leading radiation effects is more compli- 
cated than for s-channel-like processes. Finally, the run- 
ning of the electromagnetic coupling a(q^) has to be eval- 
uated in the region of small momentum transfer (q^ < 0) 
where the fit [47] of this quantity to the hadronic vacuum 
polarisation should be used. 

The Monte Carlo generator KoralW [48] has recently 
been updated to include the ISR-FSR interference effects 
as well as the proper running of a{ql). Therefore, this 
program now has reached a level of accuracy similar to the 
other state-of-the-art programs for single-W production: 
Grc4f [49], Nextcalibur [50], Swap [51], Wphact 
[52], and Wto [53]. More detailed descriptions of these 
codes can be found in Ref. [38]. It should be kept in mind 
that none of these calculations includes non-universal elec- 
troweak corrections, leading to a theoretical uncertainty 
of about ~ 5% in cross-section predictions. Although the 
final solution for a high-energy LC certainly requires a 
full O(o) calculation of the 4/-production process, a first 
step of improvement could be done by a careful expansion 
about the propagator poles of the photon and W boson. The 
electroweak 0(a) corrections to the process ej — > v c W, 
which are known [54], represent a basic building block in 
this calculation. 

Technical progress on radiative corrections to 
multi-particle production processes 

One-loop integrals become more and more cumbersome 
if the number N of external legs in diagrams increases. 
For TV > 4, however, not all external momenta are lin- 
early independent because of the four-dimensionality of 
space-time. As known for a long time [55], this fact opens 
the possibility to relate integrals with N > 4 to integrals 
with N < 4. In recent years, various techniques for ac- 
tual evaluations of one-loop integrals with N = 5, 6 have 
been worked out [56, 57] (see also references therein for 
older methods and results). The major complication in the 
treatment of 2 — > 3 processes at one loop concerns the 
numerical evaluation of tensor 5-point integrals; in par- 
ticular, the occurrence of inverse Gram determinants in 
the usual Passarino-Veltman reduction to scalar integrals 
leads to numerical instabilities at the phase-space bound- 
ary. A possible solution to this problem was worked out in 
Ref. [57] where the known direct reduction [55] of scalar 
5-point to 4-point integrals was generalized to tensor in- 
tegrals, thereby avoiding the occurrence of leading Gram 
determinants completely. 

In the evaluation of real corrections, such as bremsstrah- 
lung, a proper and numerically stable separation of in- 
frared (soft and collinear) divergences represents one of 



the main problems. In the phase-space slicing approach 
(see Ref. [58] and references therein) already mentioned 
above, the singular regions are excluded from the "reg- 
ular" phase-space integration by small cuts on energies, 
angles, or invariant masses. Using factorization proper- 
ties, the integration over the singular regions can be done 
in the limit of infinitesimally small cut parameters. The 
necessary fine-tuning of cut parameters is avoided in so- 
called subtraction methods (see Refs. [59, 60, 61] and ref- 
erences therein), where a specially tuned auxiliary func- 
tion is subtracted from the singular integrand in such a way 
that the resulting integral is regular. The auxiliary func- 
tion has to be chosen simple enough, so that the singular 
regions can be integrated over analytically. In Ref. [59] the 
so-called "dipole subtraction approach" has been worked 
out for massless QCD. The technique admits a convenient 
construction of such auxiliary functions for arbitrary one- 
parton emission processes, without the need of any further 
complicated analytical integrations. The dipole subtraction 
formalism was subsequently worked out for photon emis- 
sion off massive fermions in Ref. [60] and for QCD with 
massive quarks in Ref. [61]. 

Results on 2 — > 3 processes at one-loop order - 
e+e" -> vPR, ttH 

Recently some one-loop calculations of electroweak ra- 
diative corrections have been presented for 2^3 pro- 
cesses that are interesting at a future LC: c + c~ — > vDR 
[63, 64] and e+c~ -> ttH [65, 66, 67]. The results of 
Refs. [63, 66] were obtained with the GRACE-LOOP [68] 
system (see below). In Refs. [64, 65, 67] the technique [57] 
for treating tensor 5-point integrals was employed. While 
Refs. [63, 65, 66] make use of the slicing approach for 
treating soft-photon emission, the results of Refs. [64, 67] 
have been obtained by dipole subtraction and checked by 
phase-space slicing for soft and collinear bremsstrahlung. 

In e + e~ annihilation there are two main production 
mechanisms for the SM Higgs boson. In the Higgs- 
strahlung process, c + c~ — ► ZH, a virtual Z boson decays 
into a Z boson and a Higgs boson. The corresponding cross 
section rises sharply at threshold (y/s > Mz + Mr) to a 
maximum at a few tens of GeV above Mz + Mr and then 
falls off as s _1 , where y/s is the CM energy of the e + e~ 
system. In the W-boson-fusion process, e + e _ — > i/ c v c H, 
the incoming c + and c~ each emit a virtual W boson which 
fuse into a Higgs boson. The cross section of the W-boson- 
fusion process grows as In s and thus is the dominant pro- 
duction mechanism for y/s 3> Mr. At the one-loop level, 
first the contributions of fermion and sfermion loops in the 
Minimal Supersymmetric Standard Model (MSSM) have 
been evaluated in Ref. [62]. A complete calculation of the 
0{a) electroweak corrections to c + c~ — ► vvYi in the SM 
has subsequently been performed in Refs. [63, 64]. 3 Some 

3 Analytical results for the one-loop corrections have been obtained in 
Ref. [69] as MAPLE output, but a numerical evaluation of these results 
is not yet available. 



Table 1: Comparison of lowest-order cross sections for 
e + e~ — ► vvR, (<Tt ree ), of one-loop-corrected cross sec- 
tion (a), and of relative corrections (5 = er/otrcc — 1) be- 
tween Refs. [63, 64] at y/s = 500 GeV (input parameters 
of Ref. [63]) 



]M H [GeVl 


<t+ [fbl 


a ffbl 

V [1UJ 


5\%\ 


Ref. 


150 


61.074(7) 


60.99(7) 


-0.2 


[63] 




61.076(5) 


60.80(2) 


-0.44(3) 


[64] 


200 


37.294(4) 


37.16(4) 


-0.4 


[63] 




37.293(3) 


37.09(2) 


-0.56(4) 


[64] 


250 


21.135(2) 


20.63(2) 


-2.5 


[63] 




21.134(1) 


20.60(1) 


-2.53(3) 


[64] 


300 


10.758(1) 


10.30(1) 


-4.2 


[63] 




10.7552(7) 


10.282(4) 


-4.40(3) 


[64] 



Table 2: Comparison of lowest-order cross sections for 
e + c~ — > ttH, (atree), of one-loop-corrected cross section 
(a), and of relative corrections (5 = a/cr tlee — 1) between 
Refs. [66, 67] for M H = 120 GeV (input parameters of 
Ref. [66], results taken from Ref. [67]) 



Vi[GeV] 


a tree [fb} <r[fb] 5[%\ 


Ref. 


600 


1.7293(3) 1.738(2) 0.5 
1.7292(2) 1.7368(6) 0.44(3) 


[66] 
[67] 


800 


2.2724(5) 2.362(4) 3.9 
2.2723(3) 2.3599(6) 3.86(2) 


[66] 
[67] 


1000 


1.9273(5) 2.027(4) 5.2 
1.9271(3) 2.0252(5) 5.09(2) 


[66] 
[67] 



results of Refs. [63,64] are compared in Table 1. The 
agreement of the correction is within 0.2% or better w.r.t. 
the lowest-order cross sections. 

The Yukawa coupling of the top quark could be mea- 
sured at a future LC with high energy and luminosity at the 
level of ~ 5% [1] by analyzing the process c + c~ — > ttH. 
A thorough prediction for this process, thus, has to con- 
trol QCD and electroweak corrections. Table 2 summa- 
rizes some results on the electroweak 0(a) corrections of 
Refs. [66, 67]. The agreement within <~ 0.1% also holds 
for other energies and Higgs-boson masses. The results of 
the previous calculation [65] roughly agree with the ones 
of Refs. [66, 67] at intermediate values of y/s and Mh, but 
are at variance at high energies (TeV range) and close to 
threshold (large M H ). 



EVENT GENERATORS FOR 
MULTI-PARTICLE FINAL STATES 

Multi-purpose generators at parton level 

The large variety of different final states for multi- 
particle production renders multi-purpose Monte Carlo 
event generators rather important, i.e. generators that de- 
liver an event generator for a user-specified (as much as 
possible) general final state based on full lowest-order am- 
plitudes. As results, these tools yield lowest-order pre- 
dictions for observables, or more generally Monte Carlos 
samples of events, that are improved by universal radiative 
corrections, such as initial-state radiation at the leading- 
logarithmic level or beamstrahlung effects. Most of the 
multi-purpose generators are also interfaced to parton- 
shower and hadronization programs. The generality ren- 
ders these programs, however, rather complex devices and, 
at present, they are far from representing tools for high- 
precision physics, because non-universal radiative correc- 
tions are not taken into account in predictions. 

The following multi-purpose generators for multi-parton 
production, including program packages for the matrix- 
element evaluation, are available: 

• AMEGIC [70]: Helicity amplitudes are automatically 
generated by the program for the SM, the MSSM, and 
some new-physics models. Various interfaces (ISR, 
PDFs, beam spectra, ISAJET, etc.) are supported. The 
phase-space generation was successfully tested for up 
to six particles in the final state. 

• GRACE [71]: The amplitudes are delivered by a built- 
in package, which can also handle SUSY processes. 
The phase-space integration is done by BASES [72]. 
Tree-level calculations have been performed for up 
to (selected) six-fermion final states. The extension 
of the system to include one-loop corrections, the 
Grace-Loop [68] program, is under construction. 

• Madevent [73] + Madgraph [74]: The Mad- 
graph algorithm can generate tree-level matrix el- 
ements for any SM process (fully supporting parti- 
cle masses), but a practical limitation is 9,999 dia- 
grams. In addition, MADGRAPH creates MADEVENT, 
an event generator for the requested process. 

• Phegas [75] + Helac [76]: The Helac program 
delivers amplitudes for all SM processes (including all 
masses). The phase-space integration done by PHE- 
GAS has been tested for selected final states with up to 
seven particles. 

• Whizard [77] + Comphep [78] / Madgraph [74] 
/ O'MEGA [79]: Matrix elements are generated by 
an automatic interface to (older versions of) COM- 
phep, Madgraph, and (the up-to-date version of) 
O'MEGA. Phase-space generation has been tested for 
most 2^6 and some 2^8 processes; unweighed 



events are supported, and a large variety of interfaces 
(ISR, beamstrahlung, Pythia, PDFs, etc.) exists. 
The inclusion of MSSM amplitudes (O'MEGA) and 
improved phase-space generation (2 — > 6) are work in 
progress. 

All but the GRACE program make use of the multi-channel 
approach for the phase-space integration. More details can 
be found in the original references. 

Tuned comparisons of different generators, both at par- 
ton and detector level, are extremely important, but be- 
come more and more laborious owing to the large variety 
of multi-particle final states. Some progress to a facilitation 
and automization of comparisons are made by MC-tester 
project [80] and Java interfaces [81]. 

Event generators and results for e + e~ — > 6f 

Particular progress was reached in recent years in the de- 
scription of six-fermion production processes. Apart from 
the multi-purpose generators listed in the previous section, 
also dedicated Monte Carlo programs and generators have 
been developed for this class of processes: 

• SlXFAP [82]: Matrix elements are provided for all 6/ 
final states (with finite fermion masses), including all 
electroweak diagrams. The generalization to QCD di- 
agrams and the extension of the phase-space integra- 
tion for all final states is in progress. 

• EETT6F [83]: Only processes relevant for tt produc- 
tion are supported (a new version includes in the 
final state and QCD diagrams); finite fermion masses 
are possible. 

• LUSIFER [84]: All 6/ final states are possible, in- 
cluding QCD diagrams with up to four quarks; rep- 
resentative results for all these final states have been 
presented. External fermions are massless. An un- 
weighting algorithm and an interface to Pythia are 
available. 

Table 3 summarizes a brief comparison of results for some 
processes e + c~ — > 6/ relevant for tt production for mass- 
less external fermions. The results reveal good agreement 
between the various programs, where minor differences are 
presumably due to the different treatments of the bottom- 
quark Yukawa coupling, which is neglected in some cases. 

A tuned comparison of results obtained with LUSIFER 
and Whizard for a large survey of 6/ final states has been 
presented in Ref. [84]. 

"True" Monte Carlo event generators 

The event generators described above work at parton 
level (partially improved by parton showers), i.e. the final- 
state particles cannot be directly identified with particles 
in a detector. For detector simulations, these parton-level 
generators have to be interfaced with parton shower and 
hadronization programs. To facilitate this interface, the 



Table 3: Comparison of lowest-order predictions for some processes c+c — ► tt — > 6fermions at y/s = 500 GeV in 
approximation of massless fermions (input parameters and cuts of Ref. [84]) 



Cf,,n [ fbl 


AMEGIC++ 


eett6f 


LUSIFER 


Phegas 


SlXFAP 


Whizard 


u e e + e~v c bb 


5.879(8) 


5.862(6) 


5.853(7) 


5.866(9) 


5.854(3) 


5.875(3) 




5.827(4) 


5.815(5) 


5.819(5) 


5.822(7) 


5.815(2) 


5.827(3) 


vTVphb 


5.809(5) 


5.807(3) 


5.809(5) 


5.809(5) 


5.804(2) 


5.810(3) 




5.800(3) 


5.797(5) 


5.800(4) 


5.798(4) 


5.798(2) 


5.796(3) 


z/ M /i + dubb 


17.209(9) 


17.213(23) 


17.171(24) 


17.204(18) 






without QCD: 


17.097(8) 


17.106(15) 


17.095(11) 


17.107(18) 


17.096(4) 


17.103(8) 



"Les Houches accord" [85] has been designed, a set of 
Fortran common blocks for the transfer of event con- 
figurations from parton level generators to showering and 
hadronization event generators. Alternatively so-called 
"true" event generators could be used, which fully in- 
clude hadronization. The following well-known "true" MC 
generators represent general-purpose tools for investigat- 
ing not only e + c~ collisions, but also lepton-hadron and 
hadron-hadron scattering: Herwig [86], ISAJET [87], and 
Pythia [88]. The programs are supported and extended 
continuously; recent upgrades and new features relevant for 
LC physics are: 



implementation of all 2 — > 
the MSSM in lowest order; 



2 scattering processes of 



• associated Higgs production (QQ^H) in the SM and 
MSSM in lowest order (the case of charged Higgs 
bosons is not yet available in ISAJET and Pythia); 

• i?-parity-violation SUSY in HERWIG; 

• MSSM with complex parameters (cMSSM) in 
Pythia; 

• inclusion of spin correlations and matrix elements for 
3- and 4-body decays in HERWIG; 

• introduction of real corrections based on matrix el- 
ements for several e + e~ processes in HERWIG and 
Pythia (see e.g. Ref. [89]). 

Among other work in progress, the implementation of NLO 
QCD corrections in "true" generators is one of the most 
pressing issues. In particular, the matching of parton show- 
ers with matrix-element calculations at NLO has to be per- 
formed carefully; first results look very promising [90]. 

Finally, it should be mentioned that the present FOR- 
TRAN versions of HERWIG and PYTHIA will be replaced 
by C++ programs in the future [91]. 

RADIATIVE CORRECTIONS IN 
SUPERSYMMETRIC THEORIES 

In order to avoid too much overlap with the reports of the 
Higgs and SUSY groups, this section is mainly restricted 
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Figure 5: Effect of the SUSY 0(cq) corrections on the 
prediction of M w in the MSSM (taken from Ref. [93]) 



to the topics that have been presented in the Loopverein 
working group. More details and references on the subject 
can be found in Ref. [5]. 



SUSY corrections to precision observables 

The confrontation of high-precision data with theoreti- 
cal predictions is, of course, also very interesting in exten- 
sions of the SM. The one-loop corrections of the MSSM to 
muon decay and to the pseudo-observables of the Z reso- 
nance have been known for many years, but not many cor- 
rections beyond one loop exist. Recently the known uni- 
versal corrections of 0(aa s ) [92] to the p-parameter have 
been supplemented by the terms of (D(a%), 0(atc*b) and 
0{al) in Ref. [93], where a t/h = h 2 t/h /(ATr) with h t/h 
denoting the top/bottom Yukawa coupling. Figure 5 illus- 
trates the effect of the 0(a% ) corrections on the prediction 
of the W-boson mass in the MSSM. The genuine MSSM 
O(af) effects modify M w at the level of 2-3 MeV. 



Mass spectra in the MSSM 

In theories with unbroken supersymmetry the fermions 
and bosons within the same multiplet have a common mass. 
In realistic theories, such as the MSSM, SUSY is broken, 
and this statement is not valid anymore. However, the 
masses of fermions or bosons within multiplets are not all 
independent, i.e. there are non-trivial relations among mass 
parameters. Since the mass spectra of SUSY theories bear 
a lot of information on the mechanism of SUSY breaking, 
precision analyses of these spectra can serve as a window 
to grand unification. 

SUSY demands (at least) two Higgs doublets to give the 
up- and down-type fermions masses. Thus, the MSSM 
predicts the existence of two charged (H*), two neutral 
scalar (h°, H°), and one neutral pseudo-scalar (A ) Higgs 
bosons. In lowest order, the Higgs masses M H ±, Mh, 
and Mh can be calculated as functions of the A -boson 
mass Ma, the ratio of Higgs-field vacuum expectation val- 
ues, tan/3 = V2/V1, and the gauge-boson masses; in par- 
ticular, the mass Mh of the lightest Higgs boson is con- 
strained to be smaller than Mz at tree level. Beyond low- 
est order, also the remaining MSSM parameters are in- 
volved in the mass relations, and Mh can reach values up 
to about 135 GcV. The status of precision calculations of 
the neutral Higgs-boson masses has been recently reviewed 
in Ref. [94]. All available corrections 4 are implemented in 
the program FeynHiggs [96]. The predictions are based 
on full one-loop calculations and on the leading effects in 
two-loop order, i.e. the corrections of the order C(a s a t ), 
O(o% ) and 0(ot s ot\,) (see Ref. [94] and references therein). 
Recently the corrections of C(aitab) and O(a^) have be- 
come available [97]. The current theoretical uncertainty in 
the Higgs-mass predictions is about 3 GeV [94], but a fur- 
ther reduction to < 0.5 GeV should be reached to match 
the accuracy needed for a LC. In this context a proper defi- 
nition of tan j3 in higher perturbative orders is crucial, since 
different renormalization schemes (see, e.g., Ref. [98]) for 
tan (3 lead to rather different relations between tan (3 and 
physical observables such as the Higgs-boson masses. 

In the sector of charginos and neutralinos of the MSSM 

only three out of the six mass parameters are independent. 

In Refs. [99, 100] the three masses m^a of the heaV- 
^2, 3, 4 

ier neutralinos have been expressed in terms of the mass 
TO-fO of the lightest neutralino and of the masses m-+ of 
the two charginos, including the complete one-loop cor- 
rections, which depend also on the other MSSM parame- 
ters. The corrections modify the calculated masses by up to 
several GeV. The on-shell renormalization in the chargino- 
neutralino sector is described in Refs. [99, 100, 101] in de- 
tail. 

The relations among sfermion masses, together with the 



It should be mentioned that the full effective potential has been pre- 
sented in Ref. [95] at the two-loop level. However, the precise relation of 
these results with the parameters defined in the on-shell renormalization 
scheme used in MSSM parameter analyses has not been worked out so 
far. 



corresponding on-shell renormalization, are worked out in 
Refs. [101, 102]. For each generation, one of the four 
squark masses and one of the three slepton masses can be 
calculated from the other sfermion masses (and the other 
MSSM parameters entering at one loop). The one-loop cor- 
rections can amount to about 5%. 

Higgs-boson and SUSY-particle decays in the 
MSSM 

Analyses of particle decays are of great importance for 
the reconstruction of coupling structures and, thus, of in- 
teraction Lagrangians. The rich particle content of SUSY 
theories leads to a large variety of decay cascades, which 
depend in detail on the chosen scenario. A discussion 
of phenomenological implications and of radiative correc- 
tions to SUSY-particle decays can be found in Ref. [5] and 
Ref. [103], respectively. 

The decay widths of Higgs bosons in the MSSM re- 
ceived much attention in recent years, so that the predic- 
tions are well elaborate. Precise predictions can be ob- 
tained with the programs FeynHiggsDecay (based on 
Ref. [104]) and HDECAY [105]. Recently the electroweak 
0(a) corrections to the decay of the CP-odd A boson into 
sfermion pairs have been calculated in Ref. [106]. For the 
Higgs decay (f> — ► bb ((j) = h°, H°, A ), which is of partic- 
ular importance for light Higgs bosons, the resummation of 
the leading SUSY-QCD effects and the related theoretical 
uncertainty have been discussed in Ref. [107] (for previous 
work on — > bb see references therein). 

Apart from considering integrated decay rates, it is in- 
teresting to inspect distributions of decay products, which 
is important for the determination of the spin and parity of 
the decaying particle. This task requires the development 
of appropriate Monte Carlo tools. For the <p — * t + t~ de- 
cay, for instance, the Tauola program was extended to 
include r-spin correlations in Ref. [108]. 

SUSY-particle production 

The direct production of SUSY particles, if they exist, is 
among the most interesting issues at future colliders. In or- 
der to determine the properties (mass, spin, decay widths, 
couplings) of these new particles, precise measurements 
and predictions of the corresponding cross sections at the 
same level of accuracy are necessary, i.e. radiative correc- 
tions have to be taken into account. 

The electroweak radiative corrections to the production 
of sfermions, c + c~ — ► //, and charginos, e + e~ — ► x + X~, 
were worked out in Refs. [109] and [110], respectively. 
Since in these calculations the sfermions and charginos 
are assumed to be stable, the results are relevant for en- 
ergies a few decay widths above the production thresh- 
old. For a threshold scan the decay of the sfermions as 
well as non-resonant coherent background effects have to 
be included; in Ref. [Ill] off-shell effects in sfermion-pair 
production have been investigated at tree level (improved 



by universal Coulomb-like corrections). Theoretically the 
whole issue is very similar to a description of the process 
e+e~ — > WW — ■> 4/ which is discussed above. 

The SUSY multiplet structure, in particular, predicts 
that the strengths of the gauge-boson-fermion and gauge- 
boson-sfermion interactions, which are equal owing to 
the gauge principle, coincide with the gaugino-sfermion- 
fermion Yukawa coupling. In order to test this relation 
in SUSY QCD the processes c + c~ — * qqg, qqg, qqg, qqg 
should be studied. In Ref. [112] the QCD and SUSY-QCD 
corrections to these processes were calculated. 

More details on radiative corrections to SUSY particle 
production and decays can be found in Ref. [103] (and ref- 
erences therein). 

Renormalization of the MSSM beyond one loop 

Beyond one loop the calculation of radiative correc- 
tions within SUSY theories is highly non-trivial, because 
there is no regularization scheme that respects gauge in- 
variance and supersymmetry at the same time. For in- 
stance, conventional dimensional regularization breaks su- 
persymmetry, while dimensional reduction, which respects 
supersymmetry, is known to be not fully consistent. In 
this situation, a mathematically convincing way to per- 
form renormalization is provided by algebraic renormal- 
ization. In this framework the symmetry identities and a 
proof of renormalizability for the MSSM have been estab- 
lished in Ref. [113]. These results can serve as a basis for 
the construction of all symmetry-restoring counterterms in 
the MSSM. 

OTHER DEVELOPMENTS 

Automization of loop calculations and mainte- 
nance of computer codes 

Once the necessary techniques and theoretical subtleties 
of a perturbative calculation are settled, to carry out the 
actual calculation is an algorithmic matter. Thus, an au- 
tomization of such calculations is highly desirable, in order 
to facilitate related calculations. Various program packages 
have been presented in the literature for automatized tree- 
level, one-loop, and multi-loop calculations. A compre- 
hensive overview was, for instance, given in Ref. [1 14]; in 
the following we have to restrict ourselves to a selection of 
topics, where the emphasis is put on recent developments. 

The generation of Feynman graphs and amplitudes is a 
combinatorial problem that can be attacked with computer 
algebra. The program packages FeynArts [7] (which has 
been extended [115] for the MSSM), Qgraf [14], DI- 
ANA [11] (based on Qgraf) are specifically designed for 
this task; also the GRACE-LOOP [68] system automatically 
generates Feynman diagrams and loop amplitudes. More- 
over, the task of calculating virtual one-loop and the cor- 
responding real-emission corrections to 2 — > 2 scattering 
reactions is by now well understood. Such calculations are 



widely automated in the packages FormCalc, combined 
with LoopTools [116], and Grace-Loop [68]. 

As an illustrating example, Table 4 provides some re- 
sults on the differential cross section for c + c~ — > tt in 
lowest order as well as including electroweak 0(a) cor- 
rections. The program FA+FC was obtained from the 
output of the FeynArts and FormCalc packages and 
makes use of the LoopTools library for the numerical 
evaluation. The TOPFIT program [117, 118] was developed 
from an algebraic reduction of Feynman graphs (delivered 
from Diana) within Form; for the numerics LoopTools 
is partially employed. More detailed comparisons between 
FA+FC and TOPFIT, including also other fermion flavours, 
can be found in Ref s. [117, 119]. The Grace-Loop result 
is completely independent of the two others. The agree- 
ment between these results reflects the enormous progress 
achieved in recent years in the automization of one-loop 
calculations. 

The Grace-Loop system has recently been used in the 
calculation of the electroweak corrections to the 2 — > 3 
processes e + e~ — > uuH, ttH [63, 66], which are discussed 
above. 

Clearly the calculation of radiative corrections is a very 
laborious task, leading to rather complex and lengthy com- 
puter codes, which should be carefully documented. The 
SANC project [120] (former CalcPHEP [121]) aims at 
providing theoretical support of this kind for future accel- 
erator experiments, using the principle of knowledge stor- 
ing. This approach is rather different from the strategy of 
automization described above, which aims at generating 
completely new programs. The Sanc program contains 
another independent calculation of the 0(a) corrections 
to c+c~ — > tt, the results of which are also included in 
Table 4. 

Numerical approaches to loop calculations 

Most of the various techniques of performing loop calcu- 
lations share the common feature that the integration over 
the loop momenta is performed analytically. This proce- 
dure leads to complications at one loop if five or more ex- 
ternal legs are involved, since both speed and stability of 
programs become more and more jeopardized. At the two- 
loop level, already the evaluation of self-energy and vertex 
corrections can lead to extremely complicated higher tran- 
scendental functions that are hard to evaluate numerically. 

An idea to avoid these complications is provided by a 
more or less purely numerical evaluation of loop correc- 
tions. There are two main difficulties in this approach. 
Firstly, the appearing ultraviolet and infrared divergences 
have to be treated and canceled carefully. Secondly, even 
finite loop integrals require a singularity handling of the in- 
tegrand near so-called particle poles, where Feynman's ie 
prescription is used as regularization. 

In Ref. [122] a method for a purely numerical evaluation 
of loop integrals is proposed. Each integral is parametrized 
with Feynman parameters and subsequently rewritten with 



Table 4: Differential cross sections for c+c — > tt for selected scattering angles at y/s = 500 GeV; input parameters are 
defined in Ref. [1 17], the soft-photon cut parameter uj/^/s is set to 10~ 5 . 



COS0 


program 












[Pb] 

Born+virt+soft 






[Pb] 

Born+virt+rcal 


Vdcos 0) 


Vdcos 6) 


Vdcos 0) 


-0.9 


FA + FC 
Grace-Loop 
Sanc 
Topfit 


0.108839194076039 
0.108839194076 
0.10883919407522 
0.108839194076039 




-0.00205485893415 
-0.002054859 
-0.00205485893360 
-0.00205485893466 


0.13206(12) 
0.13229 


0.0 


FA+FC 
Grace-Loop 
Sanc 
Topfit 


0.225470464033559 
0.225470464033 
0.22547046403258 
0.225470464033559 




-0.04321416793299 
-0.043214168 
-0.04321416793300 
-0.04321416793192 


0.23513(14) 
0.23476 


+0.9 


FA+FC 
Grace-Loop 
Sanc 
Topfit 


0.491143715767761 
0.491143715767 
0.49114371576694 
0.491143715767761 




-0.16747885864057 
-0.16747886 
-0.16747885864510 
-0.16747885863793 


0.47709(21) 
0.47768 



partial integrations. The final expression consists of a quite 
simple part containing the singular terms and another more 
complicated looking part that can be integrated numeri- 
cally. The actual application of the method to a physical 
process is still work in progress. 

Another idea was proposed in Ref. [123] and applied to 
event-shape variables in c + e~ — > 3jets in NLO. In this 
approach virtual and real corrections are added before inte- 
gration. In their sum, no soft singularities, or more gener- 
ally singularities that cancel between virtual and real cor- 
rections, appear from the beginning. Nevertheless the prob- 
lem of a stable treatment of particle poles in loop ampli- 
tudes still remains. In Ref. [123] a solution via contour 
deformations in complex integration domains is described, 
but how this procedure can be generalized is not yet clear. 

CONCLUSIONS 

In spite of the complexity of higher-order calculations 
for high-energy elementary particle reactions, there has 
been continuous progress in the development of new tech- 
niques and in making precise predictions for physics at fu- 
ture colliders. However, to be prepared for a future e + e~ 
linear collider with high energy and luminosity, such as 
TESLA, an enormous amount of work is still ahead of us. 
Full two-loop predictions for c+c~ — > 2 fermion scatter- 
ing reactions, such as the Bhabha process, or full one-loop 
calculations for c + c~ — > 4fermions are more than techni- 
cal challenges. At this level of accuracy, field-theoretical 
issues such as renormalization, the treatment of unstable 
particles, etc., demand a higher level of understanding. Of 
course, both loop calculations as well as the descriptions 
of multi-particle production processes with Monte Carlo 
techniques require and will profit from further improving 
computing devices. 

It is certainly out of question that the list of challenges 



and interesting issues could be continued at will. The way 
to a future LC will also be highly exciting in precision 
physics. 
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